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INTRODUCTION 

In  order  to  preserve  and  protect  natural  resources  within  the  National  Park  System,  it  is 
important  to  understand  the  underlying  ecological  factors  controlling  and  maintaining  critical 
resource  attributes  of  the  system.  The  factors  governing  the  ecological  components  of  stream 
systems  are  highly  variable.  In  general,  the  structure  and  function  of  ecological  communities  in 
arid  streams  is  a  function  of  the  nature  and  strength  of  those  factors.  Inventories  of  natural 
resources  addressing  management  concerns  for  preserving  natural  systems  should  include 
analysis  of  the  status  of  these  factors.  Preservation  of  the  environmental  conditions  is  critical  to 
sustaining  the  natural  condition  of  aquatic  resources. 

In  lotic  systems  the  River  Continuum  Concept  (RCC;  Vannote  et  at.  1980),  provides  a 
complex  hypothesis,  or  paradigm,  that  can  be  used  as  a  framework  to  determine  and  integrate 
critical  environmental  factors  (i.e.,  both  physical  and  chemical),  and  the  biological  components 
comprising  the  stream  ecosystem.  The  RCC  describes  how  various  physical/chemical  factors, 
and  the  resulting  biological  community  change  as  small  headwater  streams  become  a  large  river. 
The  RCC  as  proposed  in  its  original  form  delineates  the  longitudinal  succession  as  observed 
in  streams  in  the  northern  mid-western  United  States  and  in  the  northern  Cascades,  however, 
continued  comparative  studies  have  resulted  in  a  broader  expansion  of  the  original  concepts 
framed  by  the  RCC  (Minshall  et  al.  1985). 

In  this  study  the  RCC  served  as  a  framework  to  identify  critical  ecological  variables  for 
collection  and  to  interpret  data  from  the  Virgin  River.  This  resource  inventory  of  the  Virgin  River 
ecosystem  was  designed  to  provide  Zion  National  Park  with  a  database  that  can  be  used  to 
establish  the  ecological  status  of  the  aquatic  communities  in  the  river,  serve  as  a  baseline  on 
which  to  monitor  and  assess  future  change,  and  predict  changes  in  community  attributes  due 
to  alteration  of  hydrological  regimes. 
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The  RCC  has  been  widely  used  by  aquatic  ecologists  to  link  instream  community 
structure  and  function  to  basin-wide  ecological  processes.  In  its  present  form,  the  RCC 
hypothesizes  that  as  one  moves  longitudinally  downstream,  local  environmental  and  biological 
factors  change,  affecting  the  dynamics  and  flow  of  energy  in  the  form  of  organic  carbon  and 
nutrients  through  the  system.  In  response,  the  aquatic  community  adapts  to  maximize  the 
available  energy. 

Changes  in  structure  and  function  of  the  aquatic  macroinvertebrate  communities  in 
response  to  these  environmental  and  biological  changes  has  been  well  documented  for  streams 
in  the  northern  Mid-west  and  Pacific  Northwest  (Cummins  1974).  The  composition  functional 
feeding  groups  within  the  community  is  a  reflection,  at  least  in  part,  of  the  available  food 
resources.  Aquatic  invertebrates  have  developed  various  functional  feeding  adaptations  (i.e., 
shredding,  gathering,  filtering,  grazing/scraping)  based  on  the  size  (coarse  vs.  fine),  location 
(e.g.  attached,  suspended,  deposited),  and/or  form  (i.e.,  allochthonous  vs.  autochthonous)  of 
the  available  organic  material.  Shredders,  for  example,  are  adapted  to  feed  primarily  on  coarse 
particulate  organic  material  (CPOM;  >1  mm),  primarily  allochthonous  organic  material  and 
associated  microbes.  Shredders  are  most  abundant  where  allochthonous  material  is  present  (i.e, 
small,  heavily  canopied  headwater  streams).  Shredders  reduce  CPOM  to  fine  and  very-fine 
particulate  organic  material  (FPOM;  <1mm  and  >63jim,  and  VPOM;  <63jim)  which,  because 
of  the  reduced  size,  are  more  readily  transported  downstream.  Energy  sources  in  mesic  streams 
appear  to  shift  from  CPOM  in  the  headwaters  to  FPOM  and  VPOM  as  one  moves  downstream. 
In  response,  the  functional  and  structural  attributes  of  the  community  shift.  As  availability  of 
transported  FPOM  and  VPOM  increases  downstream,  collector-gatherers  and  collector-filterers 
adapted  to  collect  FPOM  and  VPOM  from  the  substrate  and  the  water  column,  become  more 
prevalent.   Scrapers  adapted  to  feed  or  graze  directly  on  periphyton  standing  crops  occur  in 


conjunction  with  stream  conditions  that  favor  development  of  autochthonous  energy  sources. 
The  RCC  contends  that  community  structure  and  function,  energy  flow,  and  other 
important  features  related  to  organic  inputs  and  processing  within  the  stream  ecosystem  change 
from  upstream  to  downstream,  regulated  by  the  environmental  and  biological  changes  along  a 
longitudinal  component  characteristic  of  stream  ecosystems.  Knowledge  of  how  certain  critical 
ecosystem  variables  shape  stream  community  structure  and  function  along  the  stream  gradient 
allow  certain  qualified  predictions  to  be  made  if  these  variables  were  to  be  changed.  For  the 
purpose  of  this  study  selected  physical,  chemical,  and  biological  variables  along  the  stream 
gradient  were  examined  to  characterize  the  stream  ecosystem  of  the  East  and  North  Forks  of 
Virgin  River  in  and  above  Zion  National  Park.  An  effort  to  determine  the  role  of  the  hydrological 
regime  in  regulating  the  existing  biological  community  is  a  special  focus  of  this  report. 
Establishing  the  role  of  the  hydrological  flow  regime  in  maintaining  the  structure  and  function  of 
the  aquatic  communities  of  the  Virgin  River  system  will  allow  inferences  to  be  made  on  the 
effects  of  departure  from  these  existing  conditions  in  the  future. 


STUDY  AREA 

The  East  and  North  Forks  of  the  Virgin  River  are  located  in  the  Wasatch  and  Uintah 
Mountain,  and  the  Colorado  Plateau  Ecoregions  of  southwestern  Utah  in  and  above  Zion 
National  Park  (Omernik  1987).  The  North  Fork  begins  south  of  the  Maragount  Plateau  at  an 
elevation  of  2,500  m  and  flows  approximately  40  km  through  a  relatively  undisturbed  upper 
catchment  and  Zion  National  Park  before  joining  with  the  East  Fork  downstream  of  Springdale, 
Utah,  at  an  elevation  of  1 ,150  m.  The  East  Fork  is  approximately  70  km  long  and  begins  south 
of  the  Sevier  River  Basin  at  elevations  around  2,500  m.  The  upper  East  Fork  flows  south  passing 
adjacent  to  the  towns  of  Glendale,  Orderville,  and  Mount  Carmel  Junction  before  flowing  through 
Zion  National  Park.  Land  use  in  the  upper  North  Fork  catchment  (above  the  Zion  Narrows)  is 
limited  to  seasonal  recreation  and  rangeland  grazing.  Access  by  road  into  this  area  is  generally 
limited  to  a  few  dry  months  during  the  summer  and  fall.  Ownership  is  a  mixture  of  private, 
Bureau  of  Land  Management,  National  Forest,  and  National  Park  Service  lands.  Land  uses  in 
the  upper  East  Fork  catchment  (upstream  from  Mt.  Carmel  Junction)  are  more  extensive  and 
including  various  agricultural  activities,  urban  development,  and  pasture  and  rangeland  grazing 
(UWRL  1976).  Much  of  the  water  from  the  upper  East  Fork  is  diverted  for  domestic  and 
municipal  use,  and/or  irrigation.  As  a  result,  flows  in  the  East  Fork  below  Mount  Carmel 
Junction  are  severally  diminished  during  dry  periods.  The  majority  of  flow  during  this  period  in 
the  lower  East  Fork  within  Zion  National  Park  is  attributable  to  accretion  of  groundwater  from  the 
Navajo  Sandstone  formation.  In  both  lower  catchments  within  Zion  National  Park  recreational 
and  associated  activities  comprise  the  primary  land  uses,  especially  along  the  North  Fork  river 
corridor. 

The  majority  of  the  annual  precipitation  in  the  region  is  from  Pacific  fronts,  December 


through  March,  and  falls  as  snow  in  the  mountains  and  rain  at  lower  elevations  (Gregory  1950). 
The  precipitation  patterns  exhibited  by  the  North  Fork  catchment  reflect  a  hydrograph  dominated 
by  spring  snow-melt,  while  in  the  East  Fork  a  rainfall  runoff  hydrograph  predominates  (Hermes 
1991).  Localized  heavy  rainfall  events  bringing  moisture  from  the  Gulf  of  Mexico  are  prevalent 
from  July  through  September.  Like  most  unregulated  arid  streams,  the  East  and  North  Forks 
are  subject  to  frequent  and  intense  flash  floods.  Summer  high  flow  events  can  be  particularly 
dramatic,  with  discharges  increasing  by  orders  of  magnitude  (Diaz  1992).  The  intensity  and 
magnitude  of  these  sudden  flows  are  highly  erosive  and  have  the  potential  to  alter  the  physical 
characteristics  of  the  stream  channel  and  bed  materials. 


METHODS 

Ten  sample  sites  were  chosen  in  reaches  that  reflect  changes  in  environmental  factors 
(i.e.,  gradient,  stream  order,  discharge,  temperature)  along  the  natural  river  continuum  in  an 
attempt  to  identify  factors  that  may  affect  the  functional  and  structural  components  of  the  aquatic 
community  (Figure  1).  Elevations,  gradients,  and  stream  orders  were  determined  using 
1 :  100,000  scale  topographic  maps.  For  stream  order  the  smallest  permanent  tributaries  shown 
on  the  map  were  designated  as  first  order  (Strahler  1957).  Discharge  was  measured  directly  in 
the  field  at  the  time  of  sampling  using  velocity  area  methodology  (Bureau  of  Reclamation  1984). 
Velocities  were  measured  at  0.6  depth  using  a  Price  AA  current  meter.  Temperature  was  also 
measured  in  the  field  at  the  time  of  sampling  using  a  hand  held  thermometer. 

Each  site  was  sampled  for  a  suite  of  chemical  and  biological  variables  following  standard 
methods  referenced  in  Table  1 .  Dissolved  oxygen  and  pH  were  measured  in  the  field  at  the  time 
of  sampling.  In  addition,  water  and  organic  samples  were  collected.  Organic  material 
transported  in  the  water  column  was  collected  by  passing  water  through  a  series  of  sieves.  The 
volume  of  water  was  recorded  and  used  to  calculate  final  concentrations.  Chemical  analyses 
(i.e.,  Kjeldahl  nitrogen,  nitrate/nitrite,  total  phosphorus,  and  organic  carbon  content  of  various 
detrital  fractions,  CPOM,  FPOM,  VPOM)  were  performed  by  the  Soils  Testing  Laboratory  at 
Colorado  State  University. 

Benthic  macroinvertebrate  communities  were  collected  with  a  Surber  sampler  equipped 
with  a  260  |im  meshed  net.  Five  Surber  samples  were  collected  in  each  of  the  predominant 
benthic  habitats  identified  at  each  site.  Benthic  habitats,  substrates  where  invertebrate 
communities  were  collected,  were  characterized  by  collecting  a  representative  sample  of  the 
material  and  measuring  the  particle  size  distribution  characterized  using  the  grade  scale 
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developed  by  American  Geophysical  Union  (Chow  1964).  In  addition,  substrates  were 
photographed  to  evaluate  substrate  changes  due  to  high  flow  events.  Photos  were  also  used 
to  document  disturbance  and  indicate  degree  of  substrate  embeddedness  by  sand,  fines,  and 
clays  after  flow  events.  No  attempt  was  made  to  quantify  substrate  types  within  a  reach. 

Table  1 .  Chemical  and  biological  variables  and  methods. 


Variable 
Chemical 
pH 

Dissolved  oxygen 
Total  phosphorous 
Nitrate/nitrate 
Total  nitrogen 
Particulate  carbon  fractions 
Coarse  CPOM  (>1.0  mm) 
Fine  FPOM  (>63  urn  *  1.0  mm) 
Very  fine  VPOM  (<63  jim) 
Biological 

Benthic  macroinvertebrates 
Fish 


Method  (Reference) 

Electrode1 

Electrode1 

Persulfate  digestion1 

Cadmium  reduction1 

Persulfate  digestion/cadmium  reduction1 

Oxidation  diffusion2 


Surber  (260  jim  mesh)3 
Electrofishing  (2  pass  removal)4 


1APHA  Standard  Methods  (1975),  2Synder  &  Trofymow  (1984),  Wetzel  &  Ukens  (1979),  and 
4Hankin  (1986) 

Benthic  invertebrates  collected  were  identified  to  lowest  possible  taxonomic  level  using 

keys  in  Edmunds  era/.  (1976),  Flint  (1984),  Flowers  and  Hilsenhoff  (1975),  Hilsenhoff  (1984), 

Hilsenhoff  and  Bilmyer  (1973),  Merritt  and  Cummins  (1984),  Pennak  (1978),  Shefter  and  Wiggins 

(1986),  and  Wiggins  (1977).    Structural  analysis  of  the  benthic  macroinvertebrate  community 

consisted  of  determinations  of  density  (i.e.,  the  total  number  of  individuals  per  unit  area), 

taxonomic  richness  (i.e.,  the  number  of  taxa  present  at  a  site).    Functional  group  analysis 

included  only  the  class  insecta,  all  other  organisms  were  grouped  as  "other".  Functional  groups 

were  assigned  to  each  insect  taxa  based  on  feeding  strategies  listed  in  Merritt  and  Cummins 

(1984). 
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Chemical,  physical,  and  biological  measurements  and  collections  were  made,  in  general, 
on  a  biweekly  schedule  over  twelve  sampling  periods  (Table  2).  Actual  date  of  collections  at 
individual  sites  are  reported  in  Appendix  A  and  B.  On  the  East  Fork,  however,  EF187  was  not 
sampled  until  the  forth  sampling  period,  beginning  June  22nd.  Because  of  the  addition  of 
EF187,  the  sampling  at  EF176  was  decreased  to  once  a  month.  EF176,  because  of  the 
diminished  flows  was  atypical  of  the  natural  systems.  EF161  was  also  sampled  monthly,  due  to 
lack  of  accessibility  to  the  site.  On  the  North  Fork,  all  sites  were  sampled  biweekly,  with  the 
exception  of  NF14  which  was  also  relatively  inaccessible.  In  addition,  NF27  was  not  sampled 
during  the  sampling  period  starting  on  July  17,  due  to  impassable  roads  limiting  accessibility. 

Table  2.  Sample  sites  and  sampling  periods  for  collections. 


Beginning  Date  of  Bi-weeklv  Sampling  PeriodSite 
5/11     5/25     6/8       6/22     7/6       7/20     8/3       8/17     8/31     9/14     9/28 


EF192 
EF187 


*********** 

******** 

CCHyg  *****  *  *  * 

EF161  ***** 

CC-JC-7  *********** 

NF34  ********** 

NF30  ********** 

fy|P2"7  ******  *** 

NP14  ***** 

fsjpr-j  p  *********  ** 


Fish  communities  were  sampled  at  four  sites  on  the  North  Fork  (NF34,  NF30,  NF27, 
NF12)  and  East  Fork  (EF192,  EF187,  EF176,  EF161)  on  July  8-9.  Fish  collections  were  made 
in  cooperation  with  the  State  of  Utah,  Department  of  Game  and  Fish,  and  research  personnel 
from  the  Department  of  Biology,  University  of  Nevada,  Las  Vegas.  Two  successive  passes  (c, 
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and  cj  of  equal  effort  were  made  using  both  backpack  and  bank  electrofishing  gear  (Hankin 
1986).  Fish  captured  were  enumerated  by  species,  weighed,  measured,  and  released. 
Population  size  (N)  was  estimated  by 

N  =  C^/fo-Gj), 

with  variance  of  the  population  estimated  determined  by 

var[N]  =  [c,Vfc +o^)]/(c, +cj\ 
and  standard  error 

se(N)  =  var[N]. 
Confidence  intervals  (95  percent)  on  the  true  population  size  [N]  was  calculated  by 

1.96xse(N) 
at  each  reach  sampled  (Seber  and  Le  Cren  1967).  Three  or  more  catches  are  required  to  test 
the  assumption  of  constant  capture  probability  (White  er  a/.  1982).    For  this  analysis  it  was 
assumed  that  all  species  present  were  equally  liable  to  be  captured  and  had  the  same  capture 
probability  on  each  pass. 

Species  biomass  was  estimated  by  multiplying  the  average  weight  of  all  the  individuals 
of  each  species  captured  within  the  sample  reach  by  the  estimated  population  size  for  each 
species.  Total  biomass  was  determined  at  each  site  was  determined  by  the  addition  of  each 
individual  specie  biomass  found  per  site.  Goodman's  rule 

var(xy)  =  x*(var  y)  +  /(var  x)  -((var  x)(var  y)) 
was  used  to  calculate  the  variance  (Seber  and  Le  Cren  1967). 
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RESULTS  AND  DISCUSSION 

Physical/Chemical  Components 

Two  distinctly  different  gemorphic  zones  with  differing  environmental  attributes  (i.e, 
gradient,  shading,  substrates),  separated  by  narrow  canyons,  were  identified  on  both  the  East 
and  North  Forks  of  the  Virgin  River.  On  both  rivers  a  marked  geomorphic  transition  is  evident 
between  the  sites  located  in  the  upper  catchments  (EF192,  EF187,  EF176,  NF34,  NF30,  NF27) 
and  those  in  the  lower  catchments  (EF161,  EF157,  NF12;  Table  3).  On  both  the  East  and  North 
Forks  the  three  upper  catchment  sites  are  separated  from  the  lower  river  sites  by  narrow 
canyons,  the  Zion  Narrows  on  the  North  Fork  and  Paraunweap  Canyon  on  the  East  Fork. 

In  general,  gradients  on  the  North  Fork  were  higher  than  those  at  comparable  stream 
orders  on  the  East  Fork  (Table  3).  The  gradient  at  NF34  (68.9  m/km)  was  the  highest,  more 
than  double  that  at  EF192  (30.3  m/km).  Gradients  at  the  other  upper  North  Fork  sites,  NF30  and 
NF27  of  29.7  and  30.3  m/km,  however,  were  comparable  with  the  gradient  at  EF192.  Gradients 
tended  to  decrease  in  a  downstream  fashion.  Gradients  of  the  lower  sites  were  less  than  those 
found  at  the  sites  in  the  upper  catchment,  with  the  exception  of  EF176.  The  gradient  at  EF176, 
8.0  m/km  was  the  lowest  of  both  rivers.  Gradients  at  the  lower  sites,  NF12,  EF161  and  EF157, 
were  relatively  similar,  9.5,  10.2  and  8.5  m/km,  respectively. 

Canopy  cover  and  shading  of  the  river  by  riparian  vegetation  decreased  downstream  as 
the  river  widened,  fitting  the  general  tenets  of  the  RCC.  Riparian  vegetation  along  the  lower  sites 
was  in  many  places  discontinuous,  often  consisting  of  a  single  line  of  trees.  Sources  of  litter 
inputs  and  shading,  as  a  result,  were  reduced.  However,  due  to  canyon  walls  shading  at  EF161 , 
NF14,  and  NF12  was  notable. 
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Diversity  of  bed  materials  increased  in  complexity  as  one  moved  downstream.  Substrates 
at  the  second  order  sites,  NF34  and  EF192,  were  dominated  by  clean,  coarse  gravel,  while 
substrates  at  the  third  order  sites,  NF30,  NF27,  EF187,  and  EF176  were  more  complex,  ranging 
from  coarse  gravel  to  cobble,  small  boulder  mixtures.  Only  at  EF187  were  substrates  embedded 
with  fine  sand  and  silt  consistently  throughout  the  study  period.  Substrates  at  the  fourth  order 
site  NF12  were  comprised  primarily  of  a  mixture  of  large  cobble  and  small  boulders.  Small  areas 
of  mixed  gravel  and  sand  were  also  found  along  the  stream  margin.  Substrates  at  the  lower  East 
Fork  sites  (EF161  and  EF157)  were  more  diverse,  again  primarily  comprised  of  cobble,  small 
boulders,  but  large  areas  of  sand  and  gravel/sand  mixtures  were  also  present.  Substrates  at 
the  lower  sites  on  both  rivers  were  commonly  embedded  by  sand  transported  by  high  flows. 
After  sand  was  deposited  it  was  gradually  transported  downstream  and  underlying  substrates 
were  again  exposed.  Dense  mats  of  filamentous  algae  developed  on  the  large  boulders  as  they 
reemerged. 

Discharges  measured  at  the  sites  during  the  study  period  are  depicted  in  Figure  2  (see 
Appendix  A  for  data).  Only  at  the  lower  North  Fork  site,  NF12,  was  a  snow-melt  hydrograph 
apparent  from  the  discharge  measurements.  On  both  forks,  discharge  increased  downstream 
due  to  the  added  flow  from  tributaries  and  groundwater  accretion.  There  was,  however,  a 
reduction  of  flow  at  EF176  due  to  withdrawal  of  water  from  the  East  Fork  for  municipal  and 
irrigation  use.  Despite  these  withdrawls,  measurements  in  the  lower  East  Fork  at  EF161  and 
EF157  showed  substantial  increased  in  flow  of  several  times  that  of  the  discharge  at  EF176. 

Because  of  the  time  interval  between  discharge  measurements,  high  flows  resulting  from 
summer  rainfall  runoff  events  were  not  evident.  Discharge  measurements  taken  every  two  weeks 
indicate  that  there  was  relatively  little  apparent  fluctuation  in  discharge  throughout  the  study 
period.  Only  a  slight  increase  in  discharge  is  apparent  in  the  lower  North  Fork  during  the  August 
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Figure  2.   Discharge  (m3/sec)  measured  at  sample  sites  on  the  East  and  North  Forks. 
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17m  sampling  period.  Hydrographs  from  the  North  Fork  near  Springdale,  Utah  (0940550),  the 
East  Fork  near  Glendale  (USGS  gage  09404450),  and  the  Virgin  River  below  the  confluence  of 
the  East  and  North  Forks  at  Virgin,  Utah  (USGS  gage  09406000),  however,  reflect  the  intensity 
of  the  flows  associated  with  the  summer  rainfall  runoff  events  during  the  sampling  period  (Figure 
3;  USGS  1986).  The  hydrographs  from  these  gages  indicate  numerous  high  flow  events  did 
occur  during  the  sample  period.  Rainfall  runoff  events  resulting  in  flows  several  times  higher  than 
baseflow  were  prevalent  in  July,  August,  and  September. 


12r 


09406000 


09405500 


09404450 


JUNE  JULY 


AUG 


SEPT 


Figure  3.   Hydrographs  from  May  to  October  1987  (Water  Year  1986)  for  the  East  Fork  near 
Glendale  (09404450),  the  North  Fork  of  the  virgin  River  near  Springdale,  Utah  (09406000)  and 
the  Virgin  River  below  the  confluence  of  the  East  and  North  Forks  at  Virgin,  Utah  (0940550). 

Substrate  disturbance  was  commonly  observed  associated  with  these  high  flows.  After 

high  flows  (indicated  in  the  field  by  above  bank  high  water  and  debris  marks),  substrates  in  the 

channel  were  disturbed  (i.e.,  algal  and  periphyton  growth  on  rock  surfaces  were  eradicated, 
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evidence  of  deposition  and  scouring  of  bed  materials  was  observed,  and  channel  morphology 
altered).  Periodic  physical  disturbance  of  substrates  was  commonly  recorded  starting  in  mid- 
July,  corresponding  to  the  high  flow  events  recorded  by  the  hydrographs  at  the  USGS  gages  on 
the  East  and  North  Forks.  The  extent  and  degree  of  disturbance  (patchiness)  varied  at  sites  and 
over  time  (Table  4).  For  example,  even  though  substrates  were  disturbed  at  NF34  and  NF12, 
no  evidence  of  disturbance  was  found  at  NF30  and  NF27  during  the  second  sample  period  in 
July.  Substrates  at  EF187  did  not  appear  disturbed  (i.e.,  remaining  embedded  with  fine  sand 
and  silt)  over  the  extent  of  the  study.  EF187  is  located  at  the  USGS  gage  site  (09404450)  on  the 
East  Fork  near  Glendale.  The  hydrograph  showed  that  very  little  variation  in  flow  occurred  (see 
Figure  3),  although  substrates  at  the  site  above  EF192  and  below  EF176  were  frequently 
disturbed  (i.e.,  dramatic  scouring  and  deposition  of  bed  material,  and  alteration  of  channel 
morphology). 

Table  4.  Reported  evidence  of  disturbance  of  substrates  due  to  high  flows.  (T)  total  or  complete 
disturbance  of  bed  materials  recorded  at  site  and  (P)  partial  or  patchy  disturbance  of  substrates 
recorded  at  site  during  the  sampling  period. 

Sampling  Periods 
6/8       6/22     7/6     7/20      8/3       8/17     8/31     9/14     9/28 

EF192  T  T 

EF187 

EF176  T 

EF161  NS*       T 

EF157  T  T 

NF34  P  P  P 

NF30  P  P 

NF27  T  NS 

NF14  NS        T 

NF12  T  T 

NS*  Site  not  sampled  during  this  period 

Water  temperatures  measured  at  the  lower  sites  on  both  forks,  in  general,  were  higher  than 
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at  the  upper  sites,  with  the  exception  of  EF176  (see  Appendix  A).  Diurnal  fluctuations  of  the 
water  temperature  at  the  lower  sites  within  the  Park,  especially  in  the  lower  East  Fork,  varied  as 
much  as  10°C  (Peterson,  1991).  Temperatures  measured  at  the  lower  East  Fork  sites  mid- 
summer were  consistently  over  20°C,  while  at  the  lower  North  Fork  site  only  once  did 
temperatures  measured  exceed  20°C.  In  addition,  temperatures  at  the  upper  sites  did  not 
exceed  20°C,  with  the  exception  of  EF176.  Temperatures  measured  at  EF176  were  regularly 
above  20°C,  and  during  mid-summer  a  high  of  27°C  was  measured.  Temperatures  above  20°C 
usually  restrict  coldwater  species,  such  as  salmonids  (Hynes  1976). 

In  general,  pH  varied  from  7.7  to  8.4  with  little  change  during  the  course  of  the  collecting 
period  among  sites,  indicating  a  well-buffered  system.  Dissolved  oxygen  varied  between  7.1  and 
9.8  mg/L  and  was  always  at  or  above  saturation  (see  Appendix  A  for  individual  chemical  values). 
These  pH  and  dissolved  oxygen  ranges  observed  do  not  pose  limitation  on  a  aquatic  biota. 

Total  nitrogen  concentrations,  comprised  of  Kjeldahl  nitrogen  and  nitrate/nitrite,  were 
higher  in  the  East  Fork  than  in  the  North  Fork  (Figure  4)  and  the  levels  at  EF176,  EF161,  and 
EF157,  downstream  of  grazing,  agriculture,  and  several  small  towns  and  communities,  were 
greater  than  the  other  East  or  North  Fork  sites.  Elevated  levels  of  nitrate/nitrite  were  also 
reported  at  the  lower  East  Fork  sites  by  Sandberg  and  Sultz  (1985).  As  most  of  the  nitrogen 
fraction  at  these  sites  is  in  the  form  of  nitrate/nitrite,  this  indicates  that  the  elevated  nitrogen 
levels  are  probably  a  result  of  upstream  land  use  along  the  stream  corridor.  Elevated  nitrogen 
concentration  at  the  other  East  and  North  Fork  sites  were  only  found  during  spring  and  mid- 
summer rainfall  runoff  events.  Fisher  (1986)  reported  that  nitrate-nitrogen  concentrations  in  flood 
waters  were  significantly  higher  than  in  free-flowing  streams  and  rivers  in  the  southwest  region 
of  the  United  States. 

Total  phosphorus  concentrations  from  the  East  and  North  Forks  showed  two  periods  with 
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elevated  concentrations,  in  May  and  again  during  mid-summer  (Figure  4).  Increases  in  total 
phosphorus  concentrations  in  the  East  and  North  Forks  appeared  to  be  related  to  high  flows 
rather  than  being  diluted  by  them.  There  was  no  evidence,  however,  of  longitudinal  increases 
in  total  phosphorus  concentrations. 

In  aquatic  systems,  nutrients  such  as  nitrogen  and  phosphorus  frequently  limit  primary 
production  (algal  growth)  and  detrital  processing.  If  N:P  ratios  are  above  12:1,  phosphorus  is 
usually  considered  to  be  the  limiting  nutrient  for  photosynthesis.  For  ratios  below  this,  nitrogen 
is  usually  limiting  (Lambou  ef  a/.  1976).  The  ratio  of  nitrogen  to  phosphorus  found  at  the  North 
Fork  and  upper  East  Fork  sites  indicated  that  nitrogen  may  be  the  limiting  primary  production, 
while  at  the  EF176,  EF161,  and  EF157  phosphorus  may  be  limiting  as  a  result  of  the  elevated 
nitrogen  levels.  We  were  unable  to  assess  the  affects  of  these  levels  on  primary  production. 
Artificial  samplers  used  to  assess  algal  standing  crops  were  repeatedly  lost  during  high  flows. 

During  the  period  of  study,  VPOM  was  the  dominant  organic  fraction  being  transported  in 
both  the  East  and  North  Forks.  Of  the  three  carbon  fractions  measured  proportionally 
concentrations  of  VPOM  were  the  largest  being  transported  in  the  North  and  East  Forks  (Figure 
5).  In  fact,  the  concentration  of  VPOM  was  much  higher  than  reported  in  many  lotic  systems 
(Hobbie  and  Likens  1973,  Naiman  1976,  1982).  The  amounts  of  FPOM  and  CPOM  measured 
were  substantially  less.  No  consistent  downstream  trends  were  noted  in  any  of  the  three  fractions 
measured.  However,  in  general,  in  the  East  Fork  the  amount  of  FPOM  measured  at  EF161  and 
EF157  was  greater  than  at  sites  upstream  throughout  the  sampling  period.  Because  of  the 
longitudinal  disruption  of  flow  at  EF176,  this  would  suggest  that  the  source  of  FPOM  is  due  to 
localized  terrestrial  inputs  and  not  upstream  processes. 

Increases  in  concentrations  of  FPOM  and  CPOM  appear  to  be  related  to  increases  in  flows 
associated  with  summer  rainfall  runoff  events.  Concentrations  of  terrestrially  derived  transported 
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organic  material  entering  the  stream  carried  by  overland  flow  are  not  reflected  by  the  maximum 
concentration  being  transported  through  the  system  due  to  time  lags  between  concentration 
peaks,  runoff  peaks,  and  collections.  Bilby  and  Likens  (1978)  found  that  maximum  FPOM 
concentrations  occurred  before  peak  discharge  and  then  dropped  rapidly.  In  addition,  they 
found  that  the  length  of  time  between  runoff  events  also  affected  concentrations  measured. 
Because  samples  were  taken  at  fixed  intervals  of  time  throughout  the  study,  it  is  doubtful  if 
collections  were  taken  the  rising  limb  of  the  storm  hydrograph.  It  is  therefore  doubtful  that 
concentrations  of  organic  matter  were  measured  when  they  were  at  their  highest  concentration. 
On  both  East  and  North  Forks  in  mid-July  at  the  onset  of  high  flow  events,  when  initial  standing 
stocks  of  terrestrial  organic  materials  stored  in  dry  channels  and  adjacent  to  the  stream  was  high 
(N.J.  Hoefs,  personal  observation),  peaks  in  organic  concentrations  measured  in  the  stream  were 
greater  than  subsequently  measured  in  August  and  September. 

The  retention  ability  of  arid  streams  is  considered  to  be  relative  low,  due  to  lack  of  physical 
retention  devices  (i.e.,  wood  debris,  backwaters,  eddies,  and  rough  bed  materials).  However, 
large  quantities  of  terrestrially  derived  FPOM  and  CPOM  was  observed  incorporated  in  sediments 
and  trapped  in  backwaters  and  eddies  immediately  following  high  flow  events  (N.J.  Hoefs, 
personal  observations).  After  the  initial  increase  in  FPOM  during  mid-summer  at  both  the  lower 
East  and  North  Fork  sites  a  relative  slow  incremental  decrease  in  FPOM  concentrations  were 
measured  throughout  August,  indicating  a  slow  release  of  retained  FPOM.  The  same  trend  was 
also  evident  of  CPOM  at  NF12.  CPOM  retention  was  much  less,  as  indicated  by  the  much 
quicker  decrease  in  concentrations  measured  over  time.  The  majority  of  the  CPOM  observed 
during  this  time  was  relatively  buoyant  and,  as  a  result,  less  readily  retained. 
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Biotic  Components 

Benthic  Macroinvertebrates  Communities 

Benthic  habitats 
Benthic  macroinvertebrate  community  collections  were  made  in  the  two  prominent  habitats 
at  each  of  the  sites,  with  the  exception  of  NF34  and  EF192  where  only  one  habitat  predominated. 
At  NF34  and  EF192  the  benthic  habitat  identified  for  collecting  macroinvertebrate  communities 
was  classified  as  coarse  gravel  (16-32mm).  Two  habitats,  both  in  "fine"  and  "coarse"  substrates 
were  identified  and  sampled  at  the  third  order  sites  on  both  the  North  Fork  (NF30  and  27)  and 
East  Fork  (EF187).  The  "fine"  habitat  was  classified  as  a  mixture  of  coarse  and  very  coarse 
gravel  (32-64mm)  and  the  "coarse"  habitat  sampled  was  dominated  primarily  by  small  cobble 
(64-250mm)  mixed  with  very  coarse  gravel  and  a  few  small  boulders  (250-500mm).  The  other 
third  order  site  on  the  East  Fork,  EF176,  was  not  typical  of  the  other  three  third  order  sites. 
Initially  only  a  "fine"  substrate  was  sampled  characterized  by  a  mixture  of  medium  (8-1 6mm)  and 
coarse  gravel,  however,  with  changing  flow  conditions  an  additional  "coarse"  sample  resembling 
the  other  third  order  sites  was  also  periodically  sampled.  In  general,  at  these  sites,  "fine"  benthic 
habitats  samples  were  predominantly  coarse  gravel,  while  "coarse"  habitats  included  a  mixture 
comprised  primarily  of  small  cobbles  interspersed  with  various  size  classes  of  gravel. 

Habitats  defined  as  "fine"  and  "coarse"  at  the  sites  at  EF161,  EF157,  NF14,  and  NF12  were 
not  comparable  to  those  in  the  upper  catchments  (EF192,  EF187,  EF176,  NF34,  NF30,  NF27). 
"Fine"  substrates  were  composed  primarily  of  mixtures  of  gravel,  however,  more  sand  (*2mm), 
and  a  few  cobbles  (64-250mm)  were  included.  The  "coarse"  habitat  samples  were  taken  from 
areas  primarily  of  cobbles  and  small  boulders  with  sand  comprising  the  majority  of  the  interstitial 
material.  The  substrates  at  the  lower  sites  became  inundated  with  sand  mid-July.  After  a  period 
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of  time,  at  the  East  Fork  sites,  the  underlying  substrates  were  re-exposed  and  samples  reflect 
recolonization  of  the  original  substrates  types.  The  substrates  at  the  North  Fork  sites,  however, 
were  never  completely  excavated.  The  samples  taken  after  mid-July  are  from  isolated  boulder 
islands  (*500mm)  surrounded  by  unstable  sands. 

Community  Structure 

Eight-six  taxa  of  aquatic  insects  and  96  total  invertebrate  taxa  were  identified  in  community 
collections  from  the  East  and  North  Forks  from  May  to  October,  1987  (Table  5).  The  aquatic 
insect  taxa  collected  were  highly  diverse  covering  the  major  orders  of  aquatic  insects.  These 
numbers  are  an  underestimation  of  the  actual  number  of  taxa  present  because  of  the  large 
number  of  taxa  that  could  not  be  identified  to  species,  primarily  in  the  dipteran  family 
Chironomidae  and  also  we  quantitatively  sampled  the  predominant  habitats  and  did  not 
exhaustively  examine  all  possible  habitats.  Taxa  richness  found  in  both  the  East  and  North 
Forks,  however,  is  comparable  to  other  arid  stream  systems  (Gray  1980). 

Taxa  richness  and  density  of  the  macroinvertebrate  communities  collected  from  the  East 
and  North  Forks  revealed  several  trends  (Figures  6-15).  On  the  East  Fork  with  the  exception 
of  EF187,  the  number  of  taxa  found  at  the  sites  rarely  exceeded  15,  whereas  in  the  upper  North 
Fork  sites  NF34,  NF30,  and  NF27,  taxa  richness  commonly  reached  20.  Numbers  found  at  the 
lower  North  Fork  sites  (NF14,  and  NF12)  were  comparable  to  numbers  found  at  the  lower  East 
Fork  sites  (EF161  and  EF157). 

In  addition,  during  the  majority  of  the  collection  periods,  taxa  richness  and  densities  of 
individuals  from  the  coarse  benthic  habitats  were  higher  than  that  found  in  fine  habitats.  Coarse 
habitats  were  comprised  of  a  greater  diversity  of  substrate  materials,  with  composites  of 
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Table  5.   Benthic  macroinvertebrate  taxa  from  the  East  and  North  Forks  of  the  Virgin  River. 


ARTHROPODA 
ARACHNIDA 
ACARINA 

HYDRACARINA 
CRUSTACEA 

AMPHIPODA 

TAUTRIDAE 

Hyaleila  azteca 
ISOPODA 

ASELUDAE 

Asellus  sp. 
OSTRACODA 
INSECTA 

COLEOPTERA 

CURCUUONIDAE 

Stenopelmus  sp. 

DRYOPIDAE 

Helichus  sp.  A 

DYTISCIDAE  L 

DYTISCIDAE  A 

Agabus  rectus  A 
Oreodytes  sp.  L 

ELMIDAE 

Cleptelmis  omata  L 
Cleptelmis  omata  A 
Heterelmis  sp.  L 
Microcylioepus  sp.  L 
Microcylloepus  sp.  A 
Narpus  sp.  L 
Optioservus  sp.  L 
Optioservus  sp.  A 
Ordobrevia  sp.  A 
Zaitzevia  parvula  L 

HYDROPHIUDAE  L 

STAPHYUNIDAE 
COLLEMBOLA 

ENTOMOBRYIDAE 
DIPTERA 

ATHERICIDAE 
Atherix  sp. 

CERATOPOGONIDAE 

CERATOPOGONIDAE  P 
Forcipomyia  sp. 

CHIRONOMIDAE  P 

DIXIDAE 

Dixa  sp. 

DOLICHOPODIDAE 

EMPIDIDAE 

EMPIDIDAE  P 

Wiedemannia  sp. 

EPHYDRIDAE 

MUSCIDAE 

Limnophora  sp. 

PSYCHODIDAE 

SARCOPHAGIDAE 

SIMUUIDAE 

Simulium  sp. 
Simulium  sp.  P 

STRATIOMYIDAE 

Caloparyphus  sp. 
Eupa/yphus  sp. 

TABANIDAE 

Chrysops  sp. 
Tabanus  sp. 


INSECTA  (cont) 

DIPTERA(cont) 

TANYDERIDAE 

Protanyderus  sp. 
TIPUUDAE 
TIPUUDAE  P 

Dicranota  sp. 

Gonomyia  sp. 

Hesperoconopa  sp. 

Hexatoma  sp. 

Ormosia  sp. 

Pedicia  sp. 

Tipula  sp. 
EPHEMEROPTERA 
BAETIDAE 

Baetis  insignificans 

Baetis  sp. 

Baetis  tricaudatus 
EPHEMERELDDAE 

Drunella  doddsi 

Dmnella  grandis 

Ephemerella  sp. 

Ephemerella  inermis 
HEPTAGENIIDAE 

Epeorus  sp. 

Heptagenia  sp. 

Rhithrogena  sp. 
LEPTOPHLEBIIDAE 

Paraleptophlebia  sp. 
TRICORYTHIDAE 

Tricorythodes  sp. 

Tricorythodes  minutus 
HEMIPTERA 

GERRIDAE 

Metrobates  sp. 
NAUCORIDAE 

Ambrysus  sp. 
VEUIDAE 

Rhagovelia  sp. 
LEPIDOPTERA 
PYRAUDAE 

Petrophila  sp. 
MEGALOPTERA 

CORYDALIDAE 

Corydalus  sp. 
ODONATA 

COENAGRIONIDAE 

Argia  sp. 
CORDULEGASTRIDAE 

Cordulegaster  sp. 
PLECOPTERA 
CAPNIIDAE 
CHLOROPERUDAE 

Suwallia  sp. 

Swettsa  sp. 
NEMOURIDAE 

Amphinemura  sp. 
PERUDAE 

Hesperoperla  pacifica 
PERLODIDAE 

Isogenoides  zionensis 

Isogenoides  elongatus 

Isoperia  sp. 
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ARTHROPODA  (cont.) 
INSECTA  (cont.) 

PLECOPTERA  (cont.) 

PTERONARCYIDAE 

Pteronarcella  sp. 
TAENIOPTERYGIDAE 

Taeniopteryx  sp. 
TRICHOPTERA 

BRACHYCENTRIDAE 

Brachycentrvs  americanus 

Micrasema  sp. 
GLOSSOSOMATIDAE 

Glossosoma  sp. 
HEUCOPSYCHIDAE 

Helicopsyche  borealis 
HYDROPSYCHIDAE 

Hydropsyche  sp. 
HYDROPTIUDAE 

Hydroptila  sp. 

Leucotrichia  sp. 

Mayatrichia  sp. 

Neotrichia  sp. 

Ochrotrichia  sp. 
LEPIDOSTOMATIDAE 

Lepidostoma  sp. 
LIMNEPHIUDAE 

Hesperophylax  magnus 

Oligophlebodes  sp. 

Psychogtypha  sp. 
PHILOPOTAMIDAE 

Wormaldia  sp. 
POLYCENTROPIDAE 

Polycentropus  sp. 
PSYCHOMYIIDAE 

Tinodes  sp. 
RHYACOPHIUDAE 

RhyacoDhila  sp. 
SERICOSTOMATIDAE 

Gumaga  griscola 

ANNELIDA 

HIRUDINEA  (6) 

RHYNCHOBDELUDA 

GLOSSIPHONIIDAE 

Helobdeila  stagnalis 
PISCICOLIDAE 

Myzobdella  lugubris 

MOLLUSCA 
GASTROPODA 

PULMONATA 

LYMNAEIDAE 
PHYSIDAE 

Physa  sp.  (2) 
PELECYPODA 

SPHAERIIDAE  (2) 

NEMATOMORPHA 


L  =  Larvae,  P  =  Pupae,  A  =  Adults 
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DATE 

Figure  6.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrate  sampled  at 

EF  192  for  each  sample  periods. 
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Figure  7.  Taxa  richness  and  density  (individuais/sq.  meter)  of  invertebrates  sampled  at 
EF  187  for  each  sample  periods.  Communities  collected  in  "fine"  benthic  habitat,  are 
indicated  by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by 
hatched  bars  (see  text  for  definition  of  habitats). 
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Figure  8.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 

EF  176  for  each  sample  periods.  Communities  collected  in  fine'  benthic  habitat,  are 

indicated  by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by  hatched 

bars  (see  text  for  definition  of  habitats).  No  invertebrates  were  found  in  August  6  samples. 
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5/20  5/29 


Figure  9.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 
EF  161  for  each  sample  periods.  Communities  collected  in  "fine"  benthic  habitat,  are 
indicated  by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by 
hatched  bars  (see  text  for  definition  of  habitats). 
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Figure  10.  Taxa  richness  and  density  (ind'rviduals/sq.  meter)  of  invertebrates  sampled  at 
EF  157  for  each  sample  periods.  Communities  collected  in  fine"  benthic  habitat,  are 
indicated  by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by  hatched 
bars  (see  text  for  definition  of  habitats).  No  invertebrates  were  found  in  August  7  samples. 
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Figure  1 1 .  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 

NF  34  for  each  sample  periods. 
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Figure  12.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 
NF  30  for  each  sample  periods.  Communities  collected  in  "fine"  benthic  habitat,  are 
indicated  by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by 
hatched  bars  (see  text  for  definition  of  habitats). 
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Figure  13.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 
NF  27  for  each  sample  periods.  Communities  collected  in  fine"  benthic  habitat,  are  indicated 
by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by  hatched  bars  (see 
text  for  definition  of  habitats).  No  invertebrates  were  found  in  August  24  samples. 
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Figure  14.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 
NF  14  for  each  sample  periods.  Communities  collected  communities  collected  in  "coarse" 
habitat  are  indicated  by  hatched  bars  (see  text  for  definition  of  habitat). 
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Figure  15.  Taxa  richness  and  density  (individuals/sq.  meter)  of  invertebrates  sampled  at 
NF12  for  each  sample  periods.  Communities  collected  in  "fine"  benthic  habitat,  are 
indicated  by  black  bars,  communities  collected  in  "coarse"  habitat  are  indicated  by 
hatched  bars.  Dotted  bars  indicate  the  predominant  habitat  type  found  after 
deposition  of  sand  material  (see  text  for  definition  of  habitats). 
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boulders,  cobbles  with  interstitial  areas  mixtures  of  various  size  classes  of  sand  and  gravel,  than 
the  fine  samples,  generally  a  homogenous  coarse  gravel  at  the  upper  sites  and  sand/gravel 
mixtures  at  the  lower  sites.  At  the  upper  sites,  however,  some  of  the  highest  measurements  of 
taxa  were  found  in  the  finer  gravel  habitats.  In  general,  both  the  highest  numbers  of  taxa  and 
densities  were  observed  during  the  first  sampling  period  in  July,  after  which  decreases  in  both 
were  associated  with  site  disturbance  associated  with  high  flow  events.  Trends  in  both  habitats 
were  generally  comparable,  however,  the  invertebrates  sampled  from  finer  substrates  appeared 
more  susceptible  to  disturbance.  As  a  result,  community  trends  reflected  in  these  habitats  in 
some  cases  were  more  dramatic. 

At  EF192  (Figure  6)  both  taxa  richness  and  densities  were  highest  in  the  sampling  period 
in  early  July.  During  the  second  sampling  period  in  July  a  partial  redistribution  and  scouring  of 
the  bed  material  (comprised  primarily  of  coarse  gravel)  was  observed.  A  decrease  in  invertebrate 
taxa  and  density  were  also  apparent,  however,  recovery  of  invertebrate  densities  to  previous 
levels  was  found  the  next  sampling  period.  A  second  high  flow  event  was  noted  during  the  last 
sampling  period  in  August,  during  which  all  the  substrates  within  the  sampling  area  were 
disturbed.  Densities  were  severely  depressed  and  did  not  recover  to  previous  levels  throughout 
the  rest  of  the  study  period.  Taxa  richness,  although  slightly  reduced  directly  after  high  flows, 
did  not  appear  to  be  as  greatly  affected  by  the  disturbance  as  density. 

Substrates  at  EF187  were  not  directly  affected  by  high  flows  over  the  study  period  (Figure 
7).  There  was  no  visible  disruption  in  the  substrate  noted,  however,  changes  in  the  degree  of 
embeddedness  over  the  study  period  were  observed.  During  high  flows  (as  indicated  by 
disturbances  at  EF192  and  EF176)  prior  to  the  second  sampling  period  in  July  the  degree  of 
embeddedness  decreased,  dropping  from  approximately  25  percent  to  less  than  5  percent. 
Invertebrate  taxa  richness  and  densities  at  this  time,  and  the  sampling  period  following  were  the 
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highest  found  at  this  site.  Density  of  invertebrates  were  twice  that  of  the  other  East  Fork  sites. 
In  August  and  September  an  increase  in  degree  of  embeddedness  again  was  observed  (reported 
exceeding  75  percent)  and  during  the  last  sampling  period  it  was  noted  that  many  of  the  gravel 
substrates  were  totally  (100%)  embedded  with  fine  materials,  primarily  clays  and  silts.  Some 
decrease  in  taxa  richness  and  substantial  decreases  in  densities  primarily  in  the  finer  substrates 
were  found. 

Invertebrates  at  EF176  were  lower  than  at  the  other  East  Fork  sites  (Figure  8).  In  August 
after  a  high  flow  event  (debris  and  water  marks  were  measured  up  to  1 .5  m  above  the  stream 
bank)  which  greatly  altered  channel  morphology  and  disturbed  the  substrates,  no  invertebrates 
were  found.  Taxa  richness  the  following  sampling  period  had  recovered  to  previous  levels  and 
by  October  densities  reached  previous  levels.  In  general,  recovery  of  taxa  richness  to  previous 
levels  was  more  rapid  than  density. 

At  the  lower  East  Fork  sites,  EF161  and  EF157,  (Figures  9  and  10)  densities  and  taxa 
richness  in  the  fine  substrate  sampled,  primarily  sand  and  gravel,  initially  were  lower  than  found 
in  the  coarse  substrate  sampled.  Sand  substrates,  due  to  its  homogeneity  and  instability, 
support  fewer  taxa  at  much  lower  densities  (Hynes  1970).  Densities  of  the  invertebrate 
communities  at  these  sites  showed  reductions  that  could  be  attributed  to  high  flows.  Taxa 
richness  and  densities  decreased  after  the  second  sampling  period  in  July  at  both  sites.  Even 
though  both  taxa  richness  and  densities  recovered  somewhat  after  high  flow  events,  taxa 
richness  at  these  sites  did  not  return  to  previous  levels.  In  addition,  recovery  was  much  quicker 
in  the  coarser  substrates  sampled. 

Trends  observed  in  invertebrate  richness  and  densities,  in  the  North  Fork  were  very  similar 
to  those  observed  in  the  East  Fork.  At  the  upper  North  Fork  sites,  NF34  and  NF27,  (Figures  1 1 
and  1 3)  reduction  in  invertebrate  densities  reflected  substrate  disturbance  associated  with  high 
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flow  events.  Taxa  richness  was  also  decreased  in  many  cases,  however,  not  as  dramatically  as 
densities.  In  contrast  to  the  East  Fork,  many  of  the  high  flow  events  only  disturbed  the  finer 
substrates  or  patches  of  substrates,  as  a  result,  the  degree  and  rates  of  recovery  varied.  The 
substrate  at  site  NF30  was  less  distrubed  due  to  the  width  of  the  canyon  which  allowed  flood 
waters  to  spread.  When  sites  were  only  partially  disturbed,  recovery  rates  were  much  faster  with 
recolonization  as  a  result  of  migration  or  drift  from  adjacent  unaffected  substrates.  However, 
reductions  in  densities  in  some  cases  were  dramatic,  for  example,  at  NF27,  reductions  of 
invertebrates  in  fine  bed  materials  exceeded  75  percent.  Brooks  and  Boulton  (1991)  determined 
that  experimental  disturbance  reduced  up  to  83%  of  the  taxa  and  97%  of  the  density  of 
individuals  in  a  stream.  Recolonization  time  was  dependent  on  the  size  and/or  patchiness  of 
the  area  disturbed.  Fisher  (1986)  estimated  that  recolonization  time  for  macroinvertebrates  in 
Sycamore  Creek  was  at  least  four  weeks  after  significant  flood  induced  disturbance. 

NF14  had  lower  densities  and  fewer  taxa  than  the  other  North  Fork  site  (Figure  14).  This 
simplification  of  the  benthic  macroinvertebrate  community  could  not  directly  be  attributed  to  high 
flow  events.  However,  this  site  was  located  in  the  narrowest  reach  of  the  river  completely 
bounded  by  rock  walls.  There  are  three  aspects  of  this  site  that  make  it  different  than  the  other 
sites.  First,  the  river  is  physically  constrained  by  the  Narrows  which  consists  of  sandstone  cliffs 
up  to  one  hundred  meters  high  which  focus  the  force  of  the  water  through  the  narrow  canyon. 
As  a  result,  a  small  increase  in  flow  could  potentially  produce  a  major  disruption  of  the  benthic 
habitat.  Secondly,  due  to  the  high  degree  of  shading  primary  production  may  be  limited,  and 
third,  this  area  of  the  river  is  subject  to  a  high  frequency  of  physical  disturbance  due  to  hikers 
through  the  Narrows. 

At  NF12  (Figure  15)  fluctuations  in  the  invertebrate  densities  and  taxa  could  be  directly 
attributed  to  alteration  of  bed  materials  associated  with  high  flows.  During  the  second  sampling 
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period  in  July  a  dramatic  change  in  substrate  and  invertebrate  densities  were  found.  The  fine 
substrate  sampled  (mixture  of  various  size  classes  of  gravel)  and  the  coarse  substrate  (primarily 
boulders)  were  replaced  by  sand  that  had  been  transported  into  the  site  during  high  flows  in  mid- 
July.  The  number  of  taxa  was  only  slightly  reduced  by  the  substrate  change,  however,  densities 
were  drastically  reduced  and  did  not  recover  until  late  September  with  the  re-emergence  of  the 
boulder  substrates. 

Functional  Feeding  Group  Analysis 

Analysis  of  the  feeding  groups  of  the  benthic  invertebrates  before  and  after  the  onset  of 
the  summer  rainfall  runoff  events  in  mid-July  provided  additional  insight  of  the  factors  affecting 
the  invertebrate  communities  in  the  East  and  North  Forks  (Figures  16-  25).  The  functional 
feeding  group  composition  at  EF192  was  comprised  primarily  of  collector-gatherers  and 
collector-filterers  (Rgure  16).  Proportion  of  filterers  and  gatherers  within  the  community  changed 
after  mid-July.  The  percentage  of  filterers  decreased  from  22  to  8  percent  and  gatherers 
increased  from  51  to  73  percent.  As  a  whole,  however,  the  percent  of  collectors  (gatherers  and 
filterers  together)  did  not  change  dramatically.  The  proportion  of  the  EF192  community 
comprised  of  shredders,  did  not  exceed  1 0  percent.  The  percentage  of  shredders  at  EF1 92,  was 
the  highest  observed  on  both  the  East  and  North  Forks.  After  mid-July  an  increase  in  shredders 
was  evident.  Increased  concentrations  in  CPOM  were  also  measured  at  EF192  after  mid-August 
associated  with  fall  riparian  leaf  abscission.  This  increase  in  available  ailochthonous  organic 
material  may  account  for  the  increased  percentage  of  shredders  observed. 

The  functional  group  composition  found  at  NF34  (Figure  21)  was  similar  to  that  found  at 
EF192.  Collectors,  primarily  gatherers  comprised  the  majority  of  the  composition:  the 
percentage  of  filterers  found  at  NF34  was  not  as  great  as  that  found  at  EF192.  Scrappers  and 
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predators  were  also  more  prevalent  at  this  site.  The  proportion  of  shredders  at  this  headwater 
site  did  not  exceed  six  percent.  Little  change  in  composition  of  the  functional  groups  occurred 
during  the  study  period  at  this  site.  Although  disturbance  (shifting  of  bed  materials  and  scouring 
of  periphyton)  was  observed  at  both  of  these  headwater  sites  and  in  some  cases  associated  with 
reduced  densities  of  invertebrates,  the  percentage  of  functional  groups  was  not  dramatically 
altered. 

Functional  group  composition  EF187  (Figure  17),  where  no  substrate  disturbance  was 
evident  during  the  study  period,  was  comparable  to  that  found  at  EF192  and  NF34.  At  EF187, 
functional  groups  remained  relatively  consistent  temporally  during  the  length  of  the  study  and 
also  spatially  between  the  predominant  benthic  habitats  analyzed. 

This  was  not  the  case  at  North  Fork  sites,  NF30  and  NF27,  where  temporal  differences 
were  evident  (Figures  22  and  23).  Before  mid-July,  in  both  benthic  habitats  sampled  gatherers 
comprised  approximately  50  percent  of  the  community.  In  the  finer  benthic  habitat,  of  coarse 
gravel,  scrapers  and  predators  comprised  a  large  proportion  the  remaining  community 
composition.  In  the  coarser  habitat,  of  cobble  and  small  boulders,  filterers,  dominated  the 
remaining  proportion  of  community.  After  mid-July,  functional  groups  collected  in  the  two  of  the 
communities  were  remarkably  similar,  gatherers  increased  and  scrapers  decreased.  Both  of  the 
benthic  habitat  types  sampled  at  these  sites  were  altered  periodically  by  high  flow  event  which 
reduced  densities  and  taxa  richness  (see  Figures  12  and  13). 

The  functional  composition  at  EF176  was  unlike  the  others  analyzed  (Figure  18).  Before 
mid-July  the  communities  collected  in  the  predominant  habitats  were  dominated  by  non-insect 
taxa.  Gatherers,  however,  were  the  dominant  insect  functional  group  followed  by  filterers 
(comprising  20  percent)  in  both  substrates.  Scrapers,  in  relatively  low  proportions,  were  only 
found  in  the  coarse  gravel.    Substrates  were  severely  altered  in  mid-July  at  this  site,  the 
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cobble/small  boulder  substrate  was  completely  inundated  by  finer  material.  In  the  coarse  gravel 
habitat  a  dramatic  shift  in  functional  feed  groups  were  found.  Filterers  were  the  predominant 
feeding  group,  comprising  62  percent  of  the  community,  followed  by  gatherers.  Scrapers  were 
completely  eliminated  from  the  community. 

At  NF14,  within  the  Narrows  collectors  were  the  primary  functional  group  found  (Figure  24). 
Gatherers  and  filters  comprised  approximately  80  percent  of  the  community  before  mid-July  and 
60  percent  after  mid-July.  This  was  the  only  case  where  a  reduction  in  collectors  was  found  after 
mid-July.  Scrapers,  surprisingly,  comprise  a  large  proportion  of  the  community  both  before  and 
after  mid-July.  Scrapers  are  adapted  to  feed  or  graze  directly  on  periphyton  standing  crops  and 
occur  in  conjunction  with  stream  conditions  that  favor  development  of  autochthonous  energy 
sources.  The  low  densities  and  number  of  invertebrate  taxa  at  this  site  are  thought  to  be  related 
to  limited  primary  production  due  to  the  nearly  constant  shading  by  the  canyon  walls.  Even 
though  scrapers  comprised  a  large  proportion  of  the  community,  actual  densities  were  low. 

Gatherers,  filterers  and  collectors,  made  up  the  predominant  portion  of  the  functional 
groups  at  the  lower  sites  EF161,  EF157,  and  NF12  (Figures  19,  20,  and  25).  Predators  and 
scrapers  in  relatively  low  proportions,  were  the  only  other  functional  groups  represented.  With 
the  exception  of  the  cobble/small  boulder  habitat  sampled  after  mid-July  at  EF161  (Figure  19), 
the  proportion  of  collectors  (gatherers  and  filters  combined)  increased  after  mid-July.  Gatherers 
select  food  based  on  particle  size  rather  than  origin.  They  consume  particles  roughly  in 
proportion  to  the  occurrence  in  the  environment.  For  example,  in  arid  stream  immediately  after 
floods,  gatherers  consume  primarily  FPOM,  but,  switch  to  diatoms  and  C/adop/7ora-derived 
detritus  as  they  become  available  (Fisher  et  al.  1982,  Gray  1980). 

No  shredders  were  found  at  these  sites  despite  the  fact  that  all  these  sites  showed  large 
increases  in  CPOM  at  the  onset  of  the  rainfall  runoff  events  in  mid-July.  Shredding  insects  are 

51 


considered  relatively  rare  in  arid  streams,  thought  primarily  due  to  the  temporal  stochasticity  of 
organic  material  inputs  and  lack  of  retention.  Retention  time  of  CPOM  does  not  appear  to  be 
long  enough  in  arid  stream  systems  for  shredders  to  complete  their  life  cycles  (Fisher  1986). 
Also  shredders  are  large  univoltine  organisms  requiring  a  full  year  to  complete  their  life  cycle. 
This  renders  them  more  vulnerable  to  disturbance. 

Fish  Communities 
Species  distributions 
The  fish  community  in  the  East  and  North  Forks  of  the  Virgin  River  is  naturally  depauperate, 
which  is  typical  of  southwestern  arid  streams  (Deacon  and  Minckley  1974).  Ten  species 
(speckled  dace,  Rhinichthys  osculus;  flannelmouth  sucker,  Catostomus  latipinnis;  desert  sucker, 
Catostomus  clarki;  mountain  sucker,  Catostomus  platyrhynchus;  Virgin  spinedace,  Lepidomeda 
mollispinis  mollispinis;  Bonneville  cutthroat  trout,  Salmo  clarki  Utah;  rainbow  trout,  Oncorhynchus 
mykiss;  brown  trout,  Salmo  tnjtta;  brook  trout,  SaJvelinus  fontinalis  and  goldfish,  Carassius 
auratus)  were  collected  during  the  July  of  1987  (Table  6).  The  speckled  dace,  flannelmouth 
sucker,  desert  sucker,  and  Virgin  spinedace  are  all  native  to  the  Virgin  River  basin;  the  latter  an 
endemic  subspecies.  It  is  unknown  whether  Bonneville  cutthroat  or  the  mountain  sucker  are 
native  to  the  Bonneville  basin  (Cross  1975),  or  introduced  by  man,  or  arrived  naturally  as  a  result 
of  headwater  capture  from  the  Bonneville  basin  (R.J.  Behnke,  Colorado  State  University,  Fort 
Collins,  CO.  personal  communication).  If  cutthroat  trout  are  native  to  the  Virgin  River  basin, 
however,  one  would  expect  to  find  a  species  native  to  the  Colorado  River  Basin  rather  than  the 
Bonneville  basin  (Behnke  1970). 

Salmonids  dominated  the  fish  communities  at  the  upper  sites  on  both  rivers  (NF34,  NF30, 
NF27,  EF192,  and  EF187;  Table  6).  At  the  upper  North  Fork  sites  cutthroat  trout  were  dominant, 
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comprising  46  percent  of  the  individuals  present.  Brook  trout  and  rainbow  trout  were  also 
present  in  relatively  high  numbers  at  NF34,  possibly  a  result  of  stocking  to  enhance  recreational 
fishing  associated  with  the  Navajo  Lake  Recreational  Area  north  of  this  site.  Cross  (1975) 
indicated  that  until  1970  rainbow  trout  were  regularly  stocked  throughout  the  Virgin  River  basin. 
On  the  East  Fork  brown  trout  was  the  dominant  species,  comprising  100  percent  of  the 
community  at  EF192,  and  84  percent  at  EF187.  No  cutthroat  trout  and  only  relatively  low 
proportions  of  rainbow  trout  were  found  in  the  upper  catchment.  No  salmonids  were  found 
below  EF187  on  the  East  Fork.  EF176  was  subject  to  extreme  daily  temperature  fluctuations. 
Daily  temperatures  routinely  exceeded  20°C  (see  Appendix  A).  Streams  with  mean  temperatures 
above  20°C  are  considered  unsuitable  for  salmonids  (Hynes  1970).  Brown  trout,  however,  were 
found  at  the  lower  North  Fork  site,  NF12,  where  water  temperature  exceeded  20°C  only  in  mid- 
summer. Canyon  shading  and/or  groundwater  accretion  appears  to  be  great  enough  to 
maintain  summer  temperature  refugia  within  the  tolerance  range  of  brown  trout. 

The  only  non-salmonid  species  (mountain  suckers,  speckled  dace  and  goldfish)  at  the 
upper  sites  were  found  at  EF187  and  EF176.  Mountain  suckers,  comprised  approximately  nine 
percent  of  the  community  at  EF187.  The  mountain  sucker  prefer  clear,  cold  water  below  21  °C 
(Sigler  and  Miller  1963).  A  mountain  sucker  was  also  found  at  EF176,  where  water  temperatures 
regularly  exceed  21  °C.  The  majority  (83  percent)  of  the  fish  community  at  EF176  was  composed 
of  speckled  dace.  Speckled  dace  are  relatively  tolerant  of  warm  water  (Sigler  and  Miller  1963). 
As  a  result,  their  distributions  did  not  appear  to  be  limited  by  the  high  temperatures  measured 
at  this  site.  Goldfish,  a  tolerant  exotic  able  to  withstand  high  temperatures  and  degraded 
conditions,  were  also  found  at  this  site.  The  goldfish  found  appeared  to  be  recent  introductions 
and  no  evidence  of  a  sustained  population  was  evident. 

The  fish  communities  at  the  lower  sites  on  the  East  and  North  Forks,  EF157  and  NF12, 
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were  composed  primarily  of  native  species,  with  the  exception  of  brown  trout  present  at  NF12. 
The  native  fish  communities  at  these  sites  were  comprised  Virgin  spinedace,  speckled  dace, 
desert  suckers  and  flannelmouth  suckers.  Virgin  spinedace  made  up  approximately  10  percent 
of  the  community  at  both  sites.  The  Virgin  spinedace  was  generally  observed  mid-water,  in  deep 
runs  and  pools  associated  with  cover.  The  upstream  distribution  limit  of  the  Virgin  spinedace 
is  thought  to  occur  on  the  North  Fork  in  Zion  National  Park  below  the  Narrows  approximately 
2-3  kilometers  above  the  1987  collection  site.  No  Virgin  spinedace  were  observed  above  the 
Temple  of  Sinawava  or  collected  above  the  Narrows.  On  the  East  Fork,  Virgin  spinedace  were 
observed  throughout  the  1987  sampling  season  above  Dennet  Canyon  on  the  East  Fork.  It  is 
presumed  that  the  falls  within  Parunuweap  mark  the  upper  boundary  of  the  Virgin  spinedace 
distribution  on  the  East  Fork  (Dr.  James  Deacon,  University  of  Nevada,  Las  Vegas,  personal 
communication).  No  Virgin  spinedace  were  found  at  the  upper  East  Fork  sites  sampled. 

Speckled  dace  comprised  63  percent  of  the  individuals  collected  and  dominated 
community  at  EF157,  while  at  NF12  they  comprised  less  than  20  percent  of  the  individuals. 
Speckled  dace  were  found  to  be  more  widely  dispersed  throughout  the  stream  reach  than  the 
Virgin  spinedace.  Sigler  and  Miller  (1963)  reported  that  speckled  dace  have  been  found  to  be 
habitat  generalists,  living  in  a  wide  variety  of  habitats.  The  presence  of  the  speckled  dace  at 
EF176  reflects  their  tolerance  for  high  temperatures  and  an  ability  to  adapt  to  a  wide  range  of 
habitats. 

At  NF12  the  fish  community  was  dominated  by  desert  suckers,  comprising  approximately 
75  percent  of  the  individuals  collected.  Relatively  few  flannelmouth  suckers  were  found  (less 
than  2  percent  of  the  community).  At  EF157,  on  the  other  hand,  the  number  of  desert  suckers 
(11  percent)  was  not  markedly  greater  than  the  number  of  flannelmouth  suckers  collected  (16 
percent).  Desert  suckers  are  morphologically  adapted  for  scraping  algae  and  diatoms  from  rock 
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surfaces,  and  have  longer  intestinal  tracks  and  shorter  gill  rakers  than  the  flannelmouth  sucker. 
The  flannelmouth,  on  the  other  hand,  is  reported  to  be  more  generalized  in  its  food  habits, 
feeding  on  a  wide  variety  of  materials.  Cross  (1975)  observed  a  predilection  in  the  desert  sucker 
for  rocky  substrates,  while  flannelmouth  suckers  appeared  more  general  in  their  habitat  selection 
being  found  predominantly  over  sand  substrates.    At  the  time  fish  communities  were  sampled 
in  July,  the  substrate  at  the  North  Fork  site  was  predominantly  large  cobble  and  small  boulders 
covered  with  a  heavy  growths  of  filamentous  algae  and  diatoms.  The  East  Fork  substrate,  at  this 
time  was  more  diverse  with  large  expanses  of  small  boulders,  cobble,  gravel,  and  sand.  Deacon 
et  a/.  (1991)  found  that  the  North  Fork  had  more  deeper  riffle  habitat,  supporting  abundant 
standing  stocks  of  periphyton,  than  did  the  East  Fork,  thus  may  provide  more  suitable  habitat 
for  desert  sucker  and  explain  its  high  relative  abundance  in  the  North  Fork.    In  addition, 
temperature  preferences  of  these  species  vary.     Deacon  et  a/.  (1987)  found  temperature 
preferenda  of  25.9°C  for  the  flannelmouth   sucker,   and   17.5°C  for  the   desert  sucker. 
Temperatures  exceeding  20°C  were  commonly  measured  at  the  East  Fork  site.  Even  though  the 
fish  communities  at  EF157  and  NF12  were  similar,  in  regards  to  species  composition,  relative 
species  abundances  appear  to  be  related  to  available  habitat. 

Fish  Biomass 

The  greatest  fish  biomass  (g/m2)  in  both  the  upper  and  lower  catchments  were  found  on 
the  North  Fork  (Table  7).  At  the  upper  sites  biomass  was  the  highest  at  NF34,  EF192,  and 
EF187,  while  biomass  at  NF30  and  NF27  were  lower.  Fish  biomass  was  low  at  EF176,  where 
only  a  few  individuals  were  found,  possible  due  to  high  seasonal  temperatures  and  lack  of 
habitat.  At  the  upper  sites,  with  the  exception  of  EF176,  fish  biomass  did  not  appear  to  be 
related  to  available  invertebrates  or  densities  at  the  time  of  sampling,  even  though  the  fish 
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communities  were  composed  exclusively  of  invertivores.  Densities  of  invertebrates  found  at 
NF30  and  NF27  (see  Figures  12  and  13)  were  much  greater  than  found  at  NF34  and  EF192  (see 
Figures  6  and  1 1).  Fish  biomass  at  NF34  and  EF192  were  greater  than  at  NF30  and  NF27  where 
high  flows  were  associated  with  reductions  in  invertebrate  densities  later  in  the  season.  The 
prevailing  fish  population  may  be  reflecting  the  greater  instability  at  these  sites  in  previous  years, 
rather  than  current  conditions. 

Table  7.  Biomass  of  fish  communities  sampled  in  the  North  and  East  Forks  of  the  Virgin  River. 

Area        Population              Total                   Biomass 
Site  Date (m2)         estimate  N*  Biomass  g" g/m2 


NF34 

7/8 

61 

25±2 

1305±487 

21.4±8.0 

NF30 

7/8 

141 

12~ 

750 

5.3 

NF27 

7/8 

127 

6±2 

508±175 

4.0±1.4 

NF12 

7/6 

1505 

780±23 

65105±8086 

43.3  ±5.4 

EF192 

7/9 

106 

53±3 

1372±646 

12.9  ±6.0 

EF187 

7/9 

190 

86±7 

3602±1211 

19.0±6.4 

EF176 

7/9 

279 

20±5 

105±69 

0.4  ±0.2 

EF157 

7/9 

277 

263±36 

2876 ±789 

10.4  ±2.8 

"Total  fish  population  ±  Standard  error. 

'Biomass  per  site  weighted  by  individual  species  weights. 

"No  fish  caught  on  second  pass;  N  based  on  one  pass. 


Biomass  was  at  the  greatest  at  NF12  due  to  the  high  density  of  desert  suckers  found. 
Biomass  at  EF157  was  statistically  significantly  less  than  found  at  NF12.  Even  though  the  fish 
biomass  was  lower  at  the  EF157  the  presence  of  a  diverse  native  community  and  the  Virgin 
spinedace  reflect  the  integrity  of  the  site.  The  Virgin  spinedace  once  considered  common  to 
abundant  throughout  the  Virgin  River  system,  however,  has  been  restricted  to  a  fraction  of  its 
original  distribution  (Valdez  era/.  1991). 
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Two  distinct  fish  communities  were  found  in  the  East  and  North  Forks.  In  the  upper 
catchments,  the  fish  community  was  dominated  by  non-native  salmonids,  while  in  the  lower 
catchment  the  fish  community  was  composed  primarily  of  native  species.  The  structure  of  both 
of  these  communities  appear  to  be  influenced  by  environmental  conditions,  primarily  the 
stochasticity  of  flows,  and  resulting  reduction  of  available  food  resources  and  habitat.  As  a 
result,  the  effects  of  biotic  interactions  on  existing  community  structure  is  minimal  or  sporadic. 
Biotic  interactions  could  become  important  determinants  of  community  structure  if  physical 
variability  is  reduced  or  species  adapted  to  existing  physical  stochasticity  are  introduced. 
Population  densities  and  biomass  of  non-native  species  appeared  curtailed  in  the  upper 
catchment  at  sites  where  food  resources  and  habitat  were  shown  to  be  affected  by  high  flows. 
The  populations  of  native  fish  comprising  the  communities  within  Zion  National  Park,  on  the 
other  hand,  having  adapted  to  stochastic  flow  events  (Minckley  1973,  Deacon  and  Minckley 
1974)  appear  relatively  stable  (Valdaz  et  al.  1991). 
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SUMMARY  AND  CONCLUSION 

The  East  and  North  Forks  hydrologically  are  highly  stochastic,  responding  to  variable 
rainfall  runoff  events  typical  of  the  region.  Over  the  period  of  study  numerous  high  flows  of 
magnitudes  able  to  disturb  stream  substrates  and  alter  channel  morphology  were  observed. 
Both  the  concentration  of  nutrients  and  organic  material  measured  in  the  East  and  North  Forks 
appeared  to  be  a  function  of  discharge  in  that  the  concentrations  of  total  nitrogen,  total 
phosphorus  and  organic  concentrations  increased  during  high  flows.  Concentrations  in 
transported  organic  material  (CPOM  and  FPOM)  appeared  to  be  primarily  due  to  the  addition 
of  terrestrially  derived  organic  matter  entering  the  system  carried  by  overland  flow. 
Proportionally,  concentrations  of  VPOM  were  the  largest  being  transported.  Nutrients  and 
organic  material  retention,  however,  were  relatively  short  in  the  streams  as  indicated  by  the  lack 
of  persistence  of  measured  elevated  levels.  The  ability  of  the  biotic  community  to  utilize  these 
resources,  especially  organic  material  was  limited  as  a  result. 

Physical  disturbance,  usually  in  the  form  of  floods,  may  be  the  primary  determinant  of 
biotic  community  structure  in  streams  (Resh  et  a/.  1988).  The  highly  variable  hydrological 
discharge,  unstable  shifting  substrates,  and  large  amounts  of  transported  inorganic  depositional 
material  were  responsible  for  reduced  taxonomic  and  functional  complexity  of  the  benthic 
macroinvertebrate  community  of  the  East  and  North  Fork  of  the  Virgin  River.  Decreases  in 
invertebrate  densities  were  commonly  associated  with  substrate  disturbance.  Invertebrate 
densities  and  taxa  richness  appeared  to  remain  relatively  constant  with  increasing  flows  until 
flows  were  sufficiently  high  to  initiate  substrate  movement,  after  which  densities  declined  in  some 
cases  to  zero.  Density  recovery  rates  to  pre-disturbance  levels  appeared  dependent  on  the 
proportion  of  area  disturbed  (patchiness).  Invertebrate  densities  were  repeatably  reduced,  but 

59 


generally  reappeared  in  numbers  approximating  former  levels.  Taxa  richness  appeared, 
however,  relatively  resilient  to  disturbance.  Gray  (1980)  suggests  that  the  life  history,  the  rapid 
life  cycles  and  rapid  growth  of  arid  stream  invertebrates  adapt  them  to  the  highly  stochastic 
hydrologic  events  characteristic  of  arid  streams  similar  to  the  East  and  North  Forks.  The 
invertebrates  found  throughout  the  East  and  North  Forks  were  primarily  collector-gatherers  and, 
in  general,  after  the  onset  of  high  flow  in  mid-July,  these  groups  increased  proportionately  to  the 
other  functional  groups.  Collector-gatherers  appeared  to  be  able  to  utilize  the  available  food 
resources  (VPOM).  Shredders  were  primarily  confined  to  the  upper  stream  sites.  Shredders  did 
not  increase  in  response  to  the  increased  CPOM  levels  measured  at  the  lower  East  and  North 
Fork  sites  within  Zion  National  Park.  Shredders  are  usually  larger  organisms  with  a  univoltine 
life  history  and,  as  a  result,  are  more  vulnerable  to,  and  less  able  to  recover  from,  frequent 
disruption. 

The  structure  of  fish  communities  appeared  to  be  related  to  temperature  and  flow  regimes. 
Coldwater  salmonids  dominated  in  the  upper  reaches,  while  native  warmwater  species 
dominated  the  fish  communities  in  the  lower  reaches  within  Zion  National  Park.  Native  fish  were 
abundant  in  the  East  and  North  Forks  within  Zion  National  Park.  Native  fish  communities  in  the 
southwestern  United  States  have  declined  with  increased  settlement  and  development.  Most  of 
these  reductions  of  native  fish  populations  have  resulted  from  loss  of  habitat,  declines  in  water 
quality,  and  the  invasion  of  non-native  species  associated  with  the  alteration  of  the  flow  regimes 
within  the  Colorado  River  system  (Hubbs  and  Deacon  1964,  Minckley  and  Deacon  1968, 
Minckley  1973,  Deacon  1979).  Currently,  populations  of  Virgin  spinedace  and  native  fish 
communities  at  the  lower  East  and  North  Fork  sites  within  Zion  National  Park  appear  to  be 
relatively  undisturbed.  However,  there  is  evidence  that  the  red  shiner  Notropis  lutrensis,  and 
mosquitofish  Gambusia  affinis  now  present  in  the  Virgin  River  basin,  can  affect  the  distribution 
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and  abundance  of  native  fish  communities  in  arid  streams  (Minckley  and  Deacon  1968).  High 
flow  variability  and  low  habitat  heterogeneity  to  which  the  native  fishes  of  the  Virgin  River  are 
uniquely  adapted,  appear  to  have  prevented  the  invasion  and  establishment  of  the  red  shiner  and 
other  non-native  species  characteristic  of  less  variable  rivers  (Deacon  1979,  Cross  1985,  Valdez 
et  al.  1991).  The  red  shiner,  an  opportunistic  species,  adapted  to  variable  habitats  east  of  the 
Rocky  Mountains  and  already  established  in  the  lower  Virgin  River,  poses  an  expecially  serious 
threat.  In  addition,  slight  temperature  reduction  could  encourage  establishment  of  trout.  Thus, 
exclusion  of  exotic  species,  shown  to  affect  the  structure  of  native  fish  communities  in  other 
southwestern  streams,  may  depend  in  part  on  maintaining  the  natural  variability  of  flow, 
temperature,  and  water  quality  in  these  sections  of  the  North  and  East  Forks  of  the  Virgin  River 
(Deacon  1979, 1980).  It  also  depends  in  part  on  preventing  introduction  of  pre-adapted  exotics. 
The  longitudinal  patterns  observed  in  the  East  and  North  Forks  differed  from  those 
hypothesized  under  the  original  tenets  of  the  RCC  (Van note  et  al.  1980).  Stream  energetics, 
nutrient  and  organic  material  dynamics,  and  biotic  community  structure  and  function  patterns 
were  primarily  controlled  by  the  stochastic  hydrologic  regime  and  differing  geomorphic  attributes 
of  the  individual  reaches  examined,  rather  than  changing  longitudinal  pattern  in  downstream 
physical  stream  variables  (i.e.,  slope  velocity,  substrates,  depth)  and  biological  processes  (i.e., 
upstream  reduction  of  organic  material,  and  alterations  of  water  chemistry  or  nutrients). 

The  RCC  did,  however,  provide  a  framework  from  which  the  critical  ecological  variables 
controlling  the  biotic  communities  in  the  East  and  North  Forks  of  the  Virgin  River  were  identified. 
A  gradual  predicable  change  in  structure  and  function  of  stream  communities  along  the 
headwater  to  large  river  gradient  was  not  apparent.  Communities  did  not  reflect  a  longitudinal 
progression  of  nutrient  additions,  and  flow  of  energy  related  to  the  input,  availability,  and 
processing  of  organic  material  by  the  biotic  component  of  the  system.    Upstream  reaches, 
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however,  did  differ  from  downstream  reaches  in  several  ways.  These  differences  were  related 
to  a  number  of  physical  characteristics  of  the  East  and  North  Fork  stream  systems,  primarily 
temperature,  habitat  stability,  localized  nutrient  and  organic  loading  along  a  specific  reach,  and 
short  retention  times  controlled  by  the  regional  geomorphology  and  hydrological  regime.  As 
a  result,  the  structure  and  function  of  the  biotic  community  reflect  life  histories  that  allow  them 
to  adapt  to  the  variability  of  this  system. 
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APPENDIX  A 
(Chemistry  Data) 


EAST  FORK  VIRGIN  RIVER 
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APPENDIX  B 
(Invertebrate  Community  Data) 
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